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Frequency and dc magnetic field dependences of dynamic susceptibility in diluted paramagnets 
LiYF 4 :Ho 3+ have been measured at liquid helium temperatures in the ac and dc magnetic fields 
parallel to the symmetry axis of a tetragonal crystal lattice. Experimental data are analyzed in the 
framework of microscopic theory of relaxation rates in the manifold of 24 electron-nuclear sublevels 
of the lowest non-Kramers doublet and the first excited singlet in the Ho 3+ ground multiplet 5 7s 
split by the crystal field of S4 symmetry. The one-phonon transition probabilities were computed 
using electron-phonon coupling constants calculated in the framework of exchange charge model 
and were checked by optical piezospectroscopic measurements. The specific features observed in 
field dependences of the in- and out-of-phase susceptibilities (humps and dips, respectively) at the 
crossings (anti-crossings) of the electron-nuclear sublevels are well reproduced by simulations when 
the phonon bottleneck effect and the cross-spin relaxation are taken into account. 

PACS numbers: 71.15.Ap, 71.15.Mb, 71.15.Rf, 71.20.Be, 75.10.Pq 



I. INTRODUCTION 

In 1968 Hcllwcge et alii found an unusual non 
monotonous behavior of the in-phase magnetic suscep- 
tibility in parallel dc and ac magnetic fields with narrow 
maxima at the dc field values corresponding to the cross- 
ings of electron-nuclear sublevels of the ground-state of 
diluted Ho 3+ ions in LaCl3 at the temperature ~1 K 
and frequencies of the ac- field in the range 10 2 — 5.10 3 
Hz. Recently, similar data were independently obtained 
on the dynamic susceptibility of the diluted paramagnct 
LiYF4iHo 3+ (0.1%) in the region of energy level cross- 
ings in Ref. Q . This latter study of the anisotropic rare- 
earth spin dynamics in the classical regime was an ex- 
tension of sub-Kelvin magnetization measurements, pre- 
viously used to reveal the role of quantum fluctuations 
of atomic magnets in the diluted LiYF4:Ho 3+ system 3 -, 
similarly to the phenomenon of resonant tunneling of 
single molecule magnets in presence of a large uniaxial 
anisotropjA In both [l[ and Q the observed peculiari- 
ties of the susceptibility dynamics, measured at relatively 
high temperature, relate to predictions of N. Bloember- 
gen and co-workers in their classical study of the cross- 
relaxation in spin systems^: " ...the susceptibility is usu- 
ally plotted at constant frequency versus applied dc field. 
There may easily occur a maximum in this plot, because 
for certain values of the external field some pairs of lev- 
els ... may become nearly equidistant" . The purpose of 
the present paper is to investigate the microscopic ori- 
gin of the susceptibility in diluted LiYF4:Ho 3+ at liquid 
helium temperatures, taking into account the different 
effects of crystal-field, electron-phonon and hypcrfine in- 
teractions, as well as cross-relaxation processes. The ex- 



tension of this many-body classical dynamics to the low- 
temperature quantum case was previously adressed in 
Refs. HQ . Nevertheless, many-spin quantum fluctuations 
in highly anisotropic systems, such as co-tunneling pro- 
cesses, have not been discussed on microscopic grounds 
up to now. Therefore, to go further in the analysis of dy- 
namical magnetic hysteresis loops measured at sub-kelvin 
temperatures in the strongly out-of-equilibrium quantum 
regime^, a clear understanding of the classical spin dy- 
namics in LiYF4:Ho 3+ is required. The quasi-resonant 
absorption of radiofrequency (5-10 MHz) ultrasound and 
ac magnetic field energy at the crossing points in the 
CaWC>4:Ho 3+ crystal was observed in Ref.Q, and some 
specific peculiarities of 19 F nuclear relaxation rates at 
the crossing points in LiYF4:Ho 3+ were found in Ref.0]. 
All these experimental findings give evidence for essential 
variations of the relaxation rates in the electron-nuclear 
subsystem, by orders of magnitude, within the vicinity 
of energy level crossing points. In the present work, new 
measurements of frequency, temperature and dc mag- 
netic field dependences of the dynamic susceptibility in 
LiYF4 single crystals containing different concentrations 
of Ho 3+ ions have been carried out. The experimental 
data arc analyzed in the framework of the microscopic 
theory of the electron-phonon interaction in a gapped ti- 
tled system. The existing data on spectral properties of 
Ho 3+ ions in LiYF^ and additional piezospectroscopic 
studies described below allowed to obtain reliable val- 
ues of electron-phonon coupling constants and rigorous 
estimates of relaxation rates in the manifold of lower 
electron- nuclear sublevels of the ground multiplet. Sim- 
ulations based on the calculated values of kinetic param- 
eters revealed remarkable differences between the com- 
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putcd and measured dependences of the dynamic sus- 
ceptibility on the external parameters (frequency, tem- 
perature, strength of the constant magnetic field). We 
found it necessary to derive a more thorough theoreti- 
cal approach accounting for the phonon bottleneck effect 
and the cross-relaxation. A self-consistent description of 
all specific features of the measured susceptibilities has 
been achieved using parameters of the phonon and mag- 
netic subsystems which have physical meaning. The pa- 
per is arranged as follows: in the first section we derive 
a general expression of the dynamic susceptibility, then 
we describe the known spectral properties of impurity 
Ho 3+ ions in LiYF4 and calculate electron-phonon cou- 
pling constants which are used to compute the relaxation 
matrix in the space of lower 24 electron-nuclear sublevels 
of the ground multiplct in the external magnetic field. 
The next section contains results of experimental stud- 
ies, and in the last section we compare measured and 
simulated in- and out-of-phase susceptibilities. 



II. THEORY OF DYNAMIC MAGNETIC 
SUSCEPTIBILITY OF A DILUTED 
PARAMAGNET 

The dynamic susceptibility x{uj) of a single paramag- 
netic ion coupled to a phonon bath in the external mag- 
netic field B = B + Bi(w,<) , (Bi(w,t) = B? exp(-£wt), 
is determined by the following expression: 



When working with eq.© and ([3]) we implicitly sup- 
pose that the equilibrium state is established in the 
bath much faster than in the electronic subsystem. Be- 
sides, we have omitted in the r.h.s. of eq.© the 
terms — 2Im^ fe ^ n mBi(w, t)pk n /fi, which do not con- 
tribute to the linear response on the weak alternating 
field Bi(w,t). The off-diagonal elements of the relax- 
ation matrix W n k = Wk—> n are the transition proba- 
bilities induced by the electron-phonon interaction. For 
the one-phonon transitions, W m ^k = w mk[ n { UJ mk) + 1] if 
the frequency of the transition to m k > and W m ^k = 
w m kn(ujkm) if Umk < where w m k is the probability 
of the spontaneous transition, and n(w m fe)is the phonon 
occupation number. The diagonal element W nn = 
— J2k Wkn determines the lifetime of the corresponding 
state n, and the coherence decay rate equals 



Ink = -^(W nn + W kk ) + r„fc (4) 

where the first term is caused by finite lifetimes, and 
T n k determines all additional contributions to the homo- 
geneous broadening of the n — ► k transition. 

Considering the energy of interaction with the time- 
dependent field -mBi(w,t) as a perturbation, we solve 
equations of motion in the linear approximation: 

p n (t) = Pan + Ap n (w)Bi exp(-iwi) (5) 



Tr(p(t)Am a ) 
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(1) 



Here Am = m (m)o; m is the operator of the ion 
magnetic moment, and (...)o defines an average value, 
corresponding to the equilibrium single ion density ma- 
trix p (H ) = e ~ Ho/kT /Tr(e- Ho/kT ) in absence of the 
time dependent field. Time evolution of the density ma- 
trix p{t) defined in the space of eigenfunctions |fc) cor- 
responding to eigenvalues Ek of the Hamiltonian H$ of 
the unperturbed electronic (or electron-nuclear) system 
can be described by the master equation for diagonal 
elements p nn = p n and by simple exponential decay 
of non-diagonal matrix elements (validity of this "secu- 
lar " ap proximation was discussed, in particular, in Rcfs. 

E3I3): 
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(n ^ k,uj nk = (E n - E k )/h) 
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(3) 



where Ap n (0) = p 0n (m nn - (m) )/kT 

Postulating that the interaction with the phonon bath 
tends to establish an equilibrium state at the instant 
value of the magnetic field, J2k W n kPok(Ho ~ m Bi) = , 
we obtain from eqs. ©© and (01 

Ap ^ n {uj) = y^(wl + W)~lWkpAmp tPP p 0p /kT (6) 

hp 
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The dynamic susceptibility (fTJ) takes the form 
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where the well known expression for the static suscep- 
tibility 
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has been used. The first line i n eg . ([5]) coincides with 
the expression presented in Ref. [l2|- 

In the case of a finite concentration of paramagnetic 
ions, we have to take into account interactions between 
the ions (dipolc-dipole, dimcr or trimer exchange, virtual 
phonon exchange) and the finite rate of a heat flow from 
the phonon reservoir to the helium bath (the phonon bot- 
tleneck effect). These interactions between paramagnetic 
ions induce energy exchange (cross-relaxation), and the 
master equation ^fy contains additional nonlinear terms^: 
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Here W^J\ m is the probability of the simultaneous tran- 
sitions tti — ► / of one ion, and p — > n of another ion. 
According to the Fermi Golden rule, this transition prob- 
ability can be written as follow a 5 ' 13 : 



W CR , = — 

r r np,lm t 



|< n, l\H 12 \p, m >\) Av S(uj pn -w lm ) (11) 



where -ffi2is the Hamiltonian of interaction between the 
ions, and Qav means a configuration averaging over the 
distribution of the paramagnetic ions in the crystal lat- 
tice. In the particular case of dipolc-dipole interactions 
between rare-earth magnetic moments m = g,//isJ( ps 
is the Bohr magneton, J is the total angular momentum, 
gj is the Lande factor), the Hamiltonian H\i has the form 
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3(JiR)(J 2 R)/i? 2 



(12) 



where R is the vector connecting the ions. The last term 
at the right-hand side of eq. (fT0|) . (...) a d, provides an 
asymptotical approach to the equilibrium distribution of 
populations in the electronic system isolated from the 
phonon bath (in adiabatic conditions) due to the cross- 
relaxation processes^. 

Taking into account the distribution of the single ion 
energies, we can write the transition probability pip in 
the following form: 



where 5^3^,5 (w) is the cross- relaxation line shape func- 
tion, S is the average energy of the interaction in fre- 
quency units (we can consider S = (<7j/zb) 2 /?LRq with 
i?o equal to a lattice constant as a scaling factor), and 
k a f3~f8 are the average values of the corresponding dimcn- 
sionless lattice factors. The Hamiltonian (fP2|) contains 
only symmetrical products of the components of the an- 
gular moments, Ji Q J2/3, and in a general case there are 
21 independent parameters kap-yS- However, this number 
can be essentially diminished due to symmetry properties 
of a concrete lattice (similarly to the elastic compliance 
tensor). In particular, in a case of S4 local symmetry of 
paramagnetic ions, we obtain 
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+ 9?2 R ki2(\J^47 + J i-J^\ 2 ) 

+ \{J7 + -J7M™ + Ji!{J l 2+-J l 2 m )\ 2 ] 

+ 9? e R kw[(J747 - + JT-Jt) + cc] 

with J ab =< a\J\b >. Because a similar expression 
may be obtained in a case of long range interactions be- 
tween the paramagnetic ions through the field of elas- 
tic lattice deformations (H& = £ pfem „ A™C£(l)C m (2) 
where C p (i) is the spherical tensor operating in the 
space of states of the ion i, and A™£ are the coupling 
constants^)wc consider seven factors k a b introduced in 
(fT4"| as the phenomenological parameters. 

If the system is not far from equilibrium, equations 
(fT0|) can be linearized. We suppose that the adiabatic 
density matrix is characterized by a single parameter, 
the adiabatic temperature. The difference 



t: 



Z^kEnApJh -T-^ffoAmBQ 
(H o ) - (H )l 



between the lattice and adiabatic temperatures is deter- 
mined from the condition that the average values of the 
electronic energy Hq — mBi obtained with the adiabatic 
density matrix and with p(t) are equal. The master equa- 
tion in the form ([2]) and the expression ([8]) for the dy- 
namic susceptibility remain valid but with the additional 
contributions to the relaxation matrix 
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Let us now take into account the finite relaxation rate of 
the phonon subsystem^ considering equations of motion 
for the phonon occupation numbers along with the (|2|)): 



dn(uj mk ) 
dt 
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n{uj m k) - n (uj mk )} (17) 
{p m [n(ui mk ) + 1] - Pkn{uj m k)) 



here no(uj mk ) = [cTcp(huj mk / kT) — 1] _1 is the equilib- 
rium occupation number of phonons at the resonance fre- 
quency uj mk > in the band with the width Auj mk , t p ^ 
is the phonon lifetime, N is the number of paramagnetic 
ions per unit volume, P m k is the density of states of res- 
onant phonons (for the low-frequency acoustic phonons 
Pmk = 3uj^ nk /2Tr 2 v 3 , where v is the average sound ve- 
locity). The stationary solution of the linearized cou- 
pled equations ((2|) and lfT7|) brings about the same re- 
sults as the solution of eq.© does (see eq.(j8j), but with 
the renormalizcd transition probabilities induced by the 
phonon bottleneck: 



2tt 2 v 3 T ph (uj 
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)Nw mn \p 0m - por; 
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(18) 

To calculate the susceptibility, we need to know, first of 
all, elements of the relaxation matrix W, i.e., the transi- 
tion probabilities for each pair of states of a paramagnetic 
ion induced by the electron-phonon interaction. At the 
next step, the parameters of the cross-relaxation rates 
and the phonon lifetimes introduced in expressions (jT4j) 
and (fT7| . respectively, may be determined from a com- 
parison of the calculated real and imaginary parts of the 
susceptibility {51) with the experimental data. 
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FIG. 1: Optical singlet-doublet transitions in the absorption 
spectrum of impurity Ho 3+ ions in LiYF4 explored in the 
piezospectroscopic measurements. 



Figffl The spectrum consists of singlets (I 1 ! and T2) and 
non-Kramers doublets ( 1^34) corresponding to irreducible 
representations T k of the S4 group. Being interested in 
low temperature magnetic properties of the system, we 
can consider only the lower part of the energy spectrum 
described by the parameterized single ion Hamiltonian 
operating in the space of 136 products \LSJJ Z > ®\II Z > 
of the electronic and nuclear spin functions (there is only 
one Holmium isotope 165 Ho with the nuclear spin 1=7/2, 
L, S, J are the electronic orbital, spin and total angular 
moments, respectively, the lowest electronic multiplct is 
5 I 8 ): 



III. SPECTRAL PROPERTIES OF ISOLATED 
Ho 3+ IONS IN LiYF 4 

The crystal lattice of LiYF4 has a space group C\ h with 
the lattice constants a = 0.5164 nm, c = 1.0741 nmJ£ 
The Ho 3+ ions substitute for Y 3+ ions in the sites with 
the point symmetry S4, coordinates of the four nearest 
fluorine ions at the distance i?i=0.2244 nm relative to 
the Ho 3+ ion in the site (000) equal (x y z), (-x -y z), 
(x -y -z), (-x y -z), where x = (t-l/2)a, y = (l/2-p)a, 
z = -qc, and p=0.2817; t=0.1645; q=0.0813. The next 
nearest fluorine ions are at the distance i?2=0.2297 nm, 
they form a deformed tetrahedron as well with x = ta, y 
= (l/2-p)a, z = (q-l/4)c). 

The energy level pattern of the electronic 4f 10 config- 
uration of the Ho 3+ ion in L1YF4 was studied in quite a 
few work o 17 ' 18 ! 19 , fragments of this pattern are shown in 



H — Ho+H e _ph 
Here 



H =H CF + H hf +H z + Hs (19) 



H CF = a 2 BlOl + ai{BlOl 
+a 6 (B° e O° e 



B\0\ 



B 4 6 Ot 



s 6 - 4 o 6 - 4 ) 



(20) 



is the crystal field Hamiltonian ( B p are the crystal 

field parameters, O p k ^ and O p ' fc ' are the real and imagi- 
nary Stevens operators, respectively; ai = a = —1/450, 
a 4 = /3 = -1/30030, a 6 = 7 = -5/(189. 143 2 ) are the 
reduced matrix elements of the Stevens operators in the 
manifold of the pure 5 Is states), H hf is the Hamiltonian 
of the magnetic hypcrfmc interaction (the electrostatic 
quadrupole interaction can be neglected), Hz = QjPb JB 
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is the electronic Zeeman energy in the magnetic field 
B; the last two terms in (fTfJ|) . the Hamiltonian of the 
electron-phonon interaction H e _ p h and the energy of in- 
teraction with random lattice strains H$, are specified 
below. As it was shown in |9|, the submillimeter EPR 
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FIG. 2: Electron-nuclear sublevels of the Ho' i+ ground dou- 
blet Is(r34) in the magnetic field B||c. The magnetic field 
strength is given in units of AB = 23.74 mT. 

spectra of Ho 3+ ions in LiYF4 are well described by the 
effective g-factor gj =1.21 (slightly different from the 
value 5/4 for the pure 5 Is multiplet), the hyperfine con- 
stant A = 0.795 GHz, and the crystal field parameters 
given in Table 1. 

It should be noted that this set of the crystal field pa- 
rameters is related to the crystallographic system of coor- 
dinates (the quantization axis z is parallel to the crystal 
c-axis). This statement is based on the earlier studies of 
the piezo-spectroscopic and nonlinear Zeeman effects in 
LiYF4:Tm 3+ 2flj2lj ]. and the anisotropic parastriction in 
the concentrated paramagnets LiTmF^^ and LiErF^. 
The calculated hyperfine structure of the ground elec- 
tronic doublet in an external magnetic field directed 
along the crystal symmetry axis is shown in FigfJ] The 
spectrum consists of the two nearly equidistant groups 
of electron-nuclear sublevels with positive and negative 
slopes, respectively. The electron-nuclear sublevels inter- 
sect at the magnetic field values B = \m! + m| AB, where 
m,m'= I z and, to first order in A, AB = A/2gjfiB- We 
can distinguish odd (C oc id, \m' — m\ = Am = 2k) and 
even (C eve n,\m' — m\ = Am =2k+l) crossing points at 
magnetic field B = (2n + l)AB and B = 2nAB, respec- 
tively, (k,n =0,1,2,3). The mixing of the ground T 3 4 dou- 
blet with the excited singlets Ti and T2 by the magnetic 
hyperfine interaction opens the gaps at the odd cross- 
ing (anti-crossing) points with the maximum values for 
Am = 2 of an order of 0.35 GHz (sec FigJ5]). However, 
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TABLE I: Crystal field parameters and electron-deformation 
coupling constants (cm -1 ) 



in addition to the avoided level crossings with Am = 2, 
gaps of comparable values (~0.28 GHz) were observed 
at the odd crossing points for the electron-nuclear sub- 
levels with Am — in the EPR spectra of the isotopi- 
cally enriched sample 7 LiYF 4 :Ho 3+ (0.1%)2. We sup- 
pose that these splittings are induced by random strains 
due to intrinsic lattice defects. As we shall see below, the 
presence of random strains is confirmed by some specific 
features of magnetic field dependences of the dynamic 
susceptibility as well, and to account for the correspond- 
ing random crystal field effects, we introduce the Hamil- 
tonian H s = a2(B\0\ + B^O^ 2 ) with the initial val- 



L Bz 2 = ±0.4752 
to describe the 



ues of parameters i?| = ^0.4525 cm 
cm -1 which has been already used in 
Am = anti-crossings observed in the EPR spectra of 
the sample isotopically enriched in Li. 

At low temperatures interactions between Ho 3+ ions 
play an important role in relaxation processes in the cou- 
pled electron-nuclear subsystemSi 3 -^. For a pair of distant 
Ho 3+ ions, the energies of lower 256 electron-nuclear sub- 
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FIG. 3: Electron-nuclear sublevels of the ground state of a 
pair of distant holmium ions in the magnetic field B||c. Ar- 
rows show half-integer crossings 



levels versus the magnetic field are shown in Fig[3J This 
figure presents a new plot of the sum of single ion ener- 
gies related to a ground ^4 doublet where the random 
"rhombic" crystal field is taken into account. It is de- 
scribed by the Hamiltonian H$ and mainly affects the 
electron-nuclear sublevels with the nuclear spin projec- 
tions m = ±1/2. The corresponding shifts (about 0.01 
cm -1 ) cause splittings of degenerate energy levels of a 
dimer which are concentrated close to the centre of grav- 
ity of the spectrum in FigJ3] (in particular, in zero mag- 
netic field, 32 states have zero energy (8 with positive, 
8 with negative, and 16 with zero slope in the magnetic 
field). It is well seen in Fig. 3 that there are many ad- 
ditional crossings of the electron-nuclear sublevels in the 
pair spectrum as compared to the single ion spectrum 
(FigJ^l), and the most important qualitative peculiarity 
is the presence of half-integer crossings (Chi) at the mag- 
netic field values B = (n + 1/2) AB for n=0,1...6. 

Widths of transitions between the electron-nuclear 
sublevels of the ground doublet and the first excited sin- 
glet measured in EPR spectra are strongly dependent on 
the holmium concentration^ at liquid helium tempera- 
tures. Assuming a Gaussian distribution, the inhomo- 
geneous broadening can be defined by the FWHM (full 
width at half maximum) of 175 MHz2. According to 
Ref.[24], frequencies of resonance transitions of the two 
non-equivalent nearest neighbor fluorine nuclei in the 
super-hyperfmc 5 F 5 - 5 I 8 spectrum equal 29.85 and 36.86 
MHz. The corresponding widths of Ho 3+ energy levels 
due to the super-hyperfine interactions may be estimated 
as 130 MHz, and the total linewidths at low concen- 
trations (~0.1%) are close to 300 MHz. The magnetic 
dipolc-dipole interactions dominate at holmium concen- 
trations larger than 0.5%, and the observed line widths 
exceed 1000 MHz in the sample containing about 1% of 



impurity Ho 3+ ions. 



IV. MODELING OF THE ELECTRON-PHONON 
INTERACTION AND SPIN-LATTICE 
RELAXATION RATES 

We shall consider interactions of impurity Ho 3+ ions 
with acoustical phonons with energies less than 40 K at 
liquid helium temperatures (energies of optical phonons 
in LIYF4 are larger than 100 K2£). The corresponding 
phonon wave lengths exceed essentially the lattice con- 
stants, and we can use the lattice elasticity theory to 
describe corresponding dynamic lattice deformations. 

In the linear approximation, the Hamiltonian of the 
electron-phonon interaction corresponding to the modu- 
lation of the crystal field by elastic waves can be written 
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HcF + Hhf + Hz, (Ja + hi) 



here the first and second terms on the right describe the 
electroii-deformational and electron-rotational interac- 
tions, respectively. Linear combinations e\(T 3 ) of the de- 
formation tensor components (e a p) which transform ac- 
cording to irreducible representations of the factor group 
C 4 h of the lattice [e(A g ) = e zz ; e(A 2 g ) = (e xx + e yy )/2; 
e(-Bg) = e xx - e yy ; e(B 2 g ) = e xy ; e\(E g ) = e xz ; e 2 (E g ) = 
e yz ) and the rotation vector $ are linear in the phonon 
annihilation (a jg ) and creation (a^) operators. In par- 
ticular, 




a/3, 1] 



2uj(q) 



1/2 



(«, 



where 



(22) 



[a/3, a] = - (e a (jq a )q 0l3 + <Tep(jq )qoa) , a = ±1 (23) 

the sum is taken over acoustical branches j ac of the lat- 
tice vibrational spectrum, N is the number of unit cells 
having the mass m, e a (j ac (\o) are the components of the 
unit polarization vector in the elastic wave with the unit 
propagation vector qo = (\/q and the frequency ujj(q). 
The components of the rotation vector 7 are given by 
the same expression (f2"2"j) where [a/3, — 1] is substituted 
for [a/3, 1} and a / /J ^ 7. Parameters of the electron- 
deformational interaction 
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)/2 



(24) 
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can be computed if the crystal field parameters for an 
impurity rare earth ion are known as explicit functions 
of coordinates X\L ta: of the host lattice ions (the vector 
Ral defines the equilibrium position of an ion A in the 
unit cell L). The corresponding calculations were carried 
out in the framework of the exchange charge model of 
the crystal field^. Parameters of the crystal field Hamil- 
tonian are represented by a sum of two terms 



-Z\(l - Op) 



XL 

2(2p+l) 
7Rxl 



<rP> 

Kl 1 



S P (R 



XL i 



O^xl^xl) (25) 



related to the electrostatic fields of point lattice ions 
with the effective charges eZ\ and "exchange" charges 
at the neighbor ions proportional to the quadratic forms 
of the overlap integrals of the 4f-elcctron and ligand 
wave functions S s =< 4/0|n".s0 >, S a =< 4f0\n"p0 >, 
S n =< 4fl\n"pl > (only outer filled n"s 2 and n"p 6 elec- 
tronic shells of ligand ions are considered) : 



S P (R) = G S S 2 (R) + G a S 2 a (R) + (2 



P(P + 1) 
12 



)G K S 2 (R) 
(26) 

Here are the numerical coefficients 2 ^, a p arc the 
shielding constants, < r p > are the moments of the 4f 
electron charge density, 0* (•&, <p) are the homogeneous 
Stevens polynomials formed from the direction cosines 
of ligand radius-vectors (the spherical coordinates of a 
ligand are R, if relative to the rare earth ion at the 
origin). Calculations were performed with a-i = 0.579, 
ct 4 = cr 6 = 27 , < r 2 > = 0.695, < r 4 > = 1.219, < r e > 
= 4.502 (atomic units)2£, the lattice sums were computed 
by the Ewald method, the dependences of the overlap 
integrals (computed with the radial wave functions from 
Refs . [29ll30| ) on the distance R (in atomic units) between 
the ions were approximated by functions So exp(— bR) 
where 50=2.2298, 0.6535, 1.3239; 6=1.2554, 0.88259, 
1.1596 for s,a and tt bonds, respectively. The values of 
the model parameters G s = G a = 5.6, G^ = 2.85 were 
obtained from fitting the calculated crystal field energies 
to the experimental data. The model results in the set 
of crystal field parameters (Table 1, column 3) which are 
well comparable with those found from the analysis of 
the optical and EPR spectra (Table 1, column 2)2. The 
coupling constants B^ x (r j ) in the Hamiltonian of the 
electron-deformational interaction presented in Table 1 
(columns 4,5,8,10,12,13,15,16) were calculated using the 
same parameters of the model. 

Now we have in hands all the data which arc necessary 
to calculate the single phonon transition probabilities in- 
troduced in section 2: 



W„ 



j=3ac 



sin OdOdup 



(27) 



E <f\J2 Bk P^ a ^+^ a p°P 



oc/3 



pk 



+i 



hu)fm[a0, -1] y^e a/ 3 7 (J 7 + 7 7 )|m > 



(n (u> m f) + 1) 



Here e a /37 is the unit antisymmetric tensor, Vj(9<p) is the 
sound velocity in the direction of the phonon wave vector 
q determined by the angular coordinates 9 and tp . The 
integrals J2j=j a Ji a Pi a ) il^ a']dQ/4:TrVj a were computed 
inRef.[3l|. 

The analysis of the magnetic field and temperature 
dependences of the relaxation rates involved the numeri- 
cal diagonalization of the matrix of the Hamiltonian Hq 
(see eq.fUH])) for the fixed values of the magnetic field 
(the matrix was constructed in the space of 136 functions 
\ 5 Is, JJ Z > ®\Hz > °f the lowest multiplet), calculations 
of transition probabilities for each pair of the electron- 
nuclear sublevels using the corresponding cigenfunctions 
of Hq, and the diagonalization of the relaxation matrix 
W. FigfJ] (a) shows the magnetic field dependences at 
the temperature 2 K of the eigenvalues of the relaxation 
matrix defined in the subspace of sublevels of the ground 
doublet. In this case transitions are mainly induced by 
the dynamic deformations of B g symmetry between the 
sublevels with the same projections of the nuclear spin. 
However, mixing of wave functions of the ground doublet 
with the wave functions of the excited singlets allows 
transitions with the nuclear spin reverting as well. At 
high magnetic fields (B >0.18 T), there are two sets of 
solutions each consisting of eight branches correspond- 
ing to allowed (fast relaxation rates exceeding 10 3 s _1 ) 
and forbidden (slow relaxation rates less than 10 2 s _1 ) 
transitions (of course, there is always a zero solution cor- 
responding to the equilibrium state of the system). At 
lower magnetic fields we see specific variations of four up- 
per branches at the Am = anti-crossings and of three 
lower branches at the Am = 2 anti-crossings pointed 
in FigO The additional narrow peaks superimposed on 
the broad peaks are caused by more narrow Am = — 2 
anti-crossings which are shifted from the Am = 2 anti- 
crossings along the magnetic field axis because the spec- 
trum is only approximately equidistant. 

Figfflb) demonstrates a remarkably different behavior 
of the considered single ion relaxation rates at the same 
temperature, 2 K, when the transitions between the first 
excited singlet and the ground doublet are taken into 
account (these transitions are induced by lattice defor- 
mations of E g symmetry). There are again two sets of 
fast and slow relaxation rates, but all rates become larger 
by one-two orders of magnitude, the dips at the Am = 
anti-crossings disappear, and six additional peaks in the 
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V. EXPRIMENTAL 



A. Piezospectroscopic measurements 




2 3 4 5 6 7 8 9 10 11 12 

Magnetic field (B/AB) 

FIG. 4: Electron-phonon relaxation rates of a single Ho 3+ ion 
in the phonon bath at the temperature 2 K. (a) - calculated 
eigenvalues of the relaxation matrix defined in the space of 
16 electron-nuclear sublevels of the ground doublet, (b) - cal- 
culated lower relaxation rates corresponding to one-phonon 
transitions between all sublevels of the ground doublet and 
the first singlet, (c) - relative energies of the electron-nuclear 
sublevels of the first excited singlet in the magnetic field B||c. 



lower branches are observed due to the anti-crossings be- 
tween the sublevels of the first excited singlet T2 which 
are shown in FigfJJc). It follows from the results of cal- 
culations of the relaxation rates at liquid helium temper- 
atures that the frequency dependence of the out-of-phasc 
susceptibility of Ho 3+ ions should exhibit a maximum at 
frequencies close to 10 5 — 10 6 Hz contrary to the ear- 
lier observations^ of this maximum at frequencies below 
10 3 Hz. This contradiction stimulated us to test the the- 
ory by direct measurements of the electron-phonon cou- 
pling constants in the piezospectroscopic experiments de- 
scribed in the following section. 



To determine parameters of the Hamiltonian (f2"Tj) cor- 
responding to the interaction of the Ho 3+ ion with lattice 
strains of _B 9 -symmetry, we measured splittings of nine 
doublets (sec FigO} induced by the uniaxial stress 
applied in the basis plane of the lattice. The sample 
LiYF4 containing 1 at.% of holmium was grown by the 
Bridgman-Stockcbargcr method and oriented with the X- 
ray diffractometer. Experiments were performed at tem- 
peratures 4.2 K and 77 K, the pressure p up to 220 MPa 
was supplied along three different directions in the (001) 
plane. Splittings of the optical lines corresponding to 
singlet-doublet transitions in the a polarized absorption 
spectra (the wave vector of the incident light was parallel 
to the c-axis) were measured directly or from the linear 
dichroism signals^. 

The uniaxial pressure in the plane (001) at the angle 
tp relative to the a axis induces the following nonzero 
components of the deformation tensor: 



f e(4) 
I e{Bl) 



SisP 

-(S n + S 12 )p/2 

-[{Sn — S12) cos2<^ + Si6 sm2<p]p 
-(25i6 cos2ip + S*66 sin2</?)p/4 



Here Sij are the components of the elastic compliance 
tensor of L1YF4 measured in Ref. [32] , in particular, Six — 
Si 2 = 18.6, Sie = 8.01, S 66 = 57.7 (10 - 6 /MPa). In the 
linear approximation, the corresponding splitting of the 
doublet with the wave functions |T34± > equals 



A(r 34 ,^ = 2 (r 34 +| J2 B k p {B l g )e{Bl)a p O k p \T^ 

pk,i=l,2 

= p(C + D cos(4(p - (p Q )) 1/2 (28) 

Using the data obtained at three different values of the 
angle ip (see Table 2), wc determined for each doublet 
the corresponding deformation potentials C and D and 
the directions ^0 of the petals of the four-fold rosettes 
(f2"8"| . Then the set of 16 electron-deformation parame- 
ters Bp(B l g ) (p k 2 2, 2 -2, 4 2, 4 -2, 6 2, 6 -2, 6 
6, 6 -6) was varied starting from the values obtained in 
the framework of the exchange charge model (Table 1, 
columns 8,10) to obtain the best estimations of the mea- 
sured splittings. Results of the fitting procedure are given 
in Table 1 (columns 7,9), and the calculated and experi- 
mental splittings are compared in Table 2. We had only 
18 independent nonlinear equations for 16 variables. It 
should be noted that the positive direction of the b-axis 
relative to the a-axis was unknown, and signs of angles 
tp were checked in the fitting procedure. Possible errors 
in the measured splittings (in particular, due to random 
internal strains and errors in the orientation of the sam- 
ple) are rather large ( ±0.15 cm _1 /100 MPa), the linear 



N 


Energy of T34 
(cm" 1 ) 


Splitting A(tp) (cm 1 


/lOOMPa) 


ip = 5° 


ip = 30.7° 


f = 45° 


1 


±8 





0.33 (0.43) 


0.67 (0.63) 


0.53 (0.53) 


2 


Is 


72 


1.12 (1.50) 


1.50 (1.98) 


1.14 (1.57) 


3 


5 I 6 


8686 


0.62 (0.37) 


0.71 (0.35) 


0.38 (0.21) 


4 


5 I 6 


8783 


1.90 (1.43) 


2.20 (1.93) 


2.00 (1.56) 


5 


% 


11240 


0.34 (0.49) 


0.44 (0.45) 


0.34 (0.30) 


6 


% 


11248 


0.43 (0.53) 


0.50 (0.78) 


0.53 (0.68) 


7 


5 I 5 


11327 


1.45 (1.33) 


1.95 (1.77) 


1.60 (1.42) 


8 


5 F 5 


15491 


0.52 (0.72) 


0.67 (0.90 


0.62 (0.85) 


9 


5 F 5 


15620 


1.01 (1.02) 


1.18 (1.33) 


0.82 (1.05) 



x'(cm 3 /mol) 



x"(cm7mol) 



TABLE II: Splittings of non-Kramers doublets in LiYF 4 :Ho 3+ 
(1%) crystals axially compressed in the basis plane. Results 
of calculations are in brackets. 



approximation for some closely spaced doublets (having 
numbers 3,5,6 in Table 2) is not valid, and the param- 
eters Bp(Bg) have not been determined unambiguously. 
However, as it is seen in Table 1, the final set of pa- 
rameters Bp(Bg) does not differ qualitatively from the 
starting values. Thus, the model used in calculations 
of the parameters of the electron-phonon interaction has 
been additionally approved. A possible error in the es- 
timations of the matrix elements of electronic operators 
in the Hamiltonian of electron-phonon interaction (f2l"j) 
which we need to calculate the transition probabilities 
between the sublevels of the ground doublet does not ex- 
ceed 24 %. 



B. Measurements of the dynamic susceptibility 

Earlier the dynamic susceptibility of diluted 
LiYi_ c Ho c F4 crystals with holmium concentrations 
c=46 %, 16.7 % and 4.5 % was studied in Ref-HHH at 
low temperatures with main attention for the spin-glass 
behavior in a random dipolar-coupled Ising magnet. 

In the present work, the dynamic susceptibility of 
highly diluted LiYi_ c Ho c F4 single crystals was measured 
with a conventional SQUID magnetometer at frequencies 
from 20 to 1200 Hz in the temperature range 1.75-4 
K. Three samples with dimensions of 4x4x10 mm 3 were 
studied in collinear ac- and dc-magnctic fields parallel to 
the c-axis of amplitude of 4 10~ 4 T and up to 0.25 T 
respectively. The holmium concentrations in these sam- 
ples (c=0.104 %; 0.157 %; 0.27 %) were determined by 
a comparison of the measured low frequency suscepti- 
bility with the calculated single ion static susceptibility. 
The in- and out-of-phase susceptibility was measured in 
a more extended dc field, frequencies, temperatures and 
concentrations, than in Ref [2j, where the results of mea- 
surements of the in- and out-of-phase susceptibilities at 
frequencies 163 and 800 Hz in the sample with holmium 
concentration c=0.04% were published. 

The results presented in FigsJS] to [7] show the non 
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FIG. 5: (Color online) Measured and simulated magnetic field 
dependences of the ac-susceptibility in LiYF4:Ho 3+ (0.157 %) 
at different temperatures ( v= 800 Hz, B||c). 



monotonous behavior with well pronounced peaks and 
dips at different crossing points of the electron-nuclear 
energy levels. In particular, the out-of-phase suscep- 
tibility measured in the vicinity of the crossing points 
for two holmium concentrations show an inversion of 
the sign of Chi peaks (dips) with respect to the back- 
ground (FigJ5]). Fig[5]prcscnts the frequency dependences 
of the in- and out-of-phase susceptibilities at three dif- 
ferent temperatures while the temperature and concen- 
tration shifts of the out-of-phase susceptibility with fre- 
quency are given in Figj9] Note that in both cases the 
dc magnetic field was taken in between two single-ion 
electron-nuclear crossing points (38.5 mT, corresponding 
to 1.62A£>). The temperature dependences of the in- 
and out-of-phase susceptibilities at different ac frequen- 
cies are shown in FigfTUl 



VI. SIMULATIONS OF THE DYNAMIC 
SUSCEPTIBILITY IN LiYF 4 : Ho 3+ AND 
DISCUSSION 



The experimental data presented in Figsl5 ffT0l are ac- 
companied by results of the corresponding simulations. 
Calculations of the susceptibilities involved the numeri- 
cal diagonalization of the single ion Hamiltonian H$ (see 
eg . (fTi?)) ) in the space of 136 electron- nuclear states of the 
lowest 5 ig multiplet for the fixed values of the external 
magnetic field, the computing of the matrix elements of 
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FIG. 6: (Color online) . Measured and simulated magnetic 
field dependences of the ac-susceptibility at the temperature 
1.75 K in LiYF4:Ho 3+ samples with different holmium con- 
centrations x (a,b - x = 0.104 %, c,d - x = 0.27 %, v = 800 
Hz, B||c). 
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FIG. 7: (Color online). Measured and simulated mag- 
netic field dependences of the ac-susceptibility in LiYF4:Ho 3+ 
(0.104 %) in the region of high-field crossings at different tem- 
peratures {v = 1200 Hz, B||c). 



electronic operators which are necessary to construct the 
relaxation matrices (fT4"|) and and at the last step 

= Rexzz(w) and x"( w ) = I m Xzz( a -') were obtained 
at different frequencies and temperatures using the ex- 
pression ([5]). 

First of all, let us consider the frequency dependences 
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FIG. 8: (Color online). Measured and simulated frequency de- 
pendences of the dynamic susceptibility in LiYF4:Ho 3+ (0.104 
%) at different temperatures (B||e, B = 38.5 mT). 
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FIG. 9: (Color online). Measured and simulated frequency 
dependences of the out-of-phase susceptibility in LiYF4:Ho 3+ 
samples at different temperatures (B||c, B = 38.5 mT). 
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FIG. 10: (Color online). Measured (symbols) and simulated 
(solid curves) temperature dependences of the dynamic sus- 
ceptibility in LiYF4:Ho 3+ (0.104 %) at different frequencies 
of the ac magnetic field (Bj|c, B = 38.5 mT). 



of the out-of-phase susceptibilities in Fig. [5] measured 
at the dc magnetic field value of 1.62 AS, which is far 
enough from the neighboring single-ion crossing points. 
The maximum of the out-of-phase susceptibility shifts 
to higher frequencies with the increasing temperature. 
These shifts give evidence for the increasing relaxation 
rates induced by the electron-phonon interaction. How- 
ever, using the one-phonon transition probabilities calcu- 
lated with the parameters of the electron-deformational 
interaction from Table 1, we found the corresponding 
maxima at the frequencies 10 4 — 10 5 Hz (see Fig |4] with 
the plotted relaxation times), which are about two orders 
of magnitude higher than the experimental values. Even 
more, Fig[5] also shows a concentration dependence of 
the maximum of the out-of-phase susceptibility, while it 
is clear that single-ion electron-phonon interaction can- 
not produce any dependence of the relaxation rate on 
the concentration of relaxing ions. Thus, to fit the simu- 
lated frequency dependences to the experimental curves, 
we had to take into account the phonon bottleneck ef- 
fect. The concentration dependence of the renormal- 
ized transition probability (|18p is determined by a fac- 
tor K m k = T p h(uj m k)N/ Aumk- Using the average sound 
velocity v = 3.10 3 m/s22, the number of paramagnetic 
ions per unit volume N = 2c/V (here V is the unit 
cell volume which contains two Y 3+ sites accessible for 



holmium ions, c is the concentration determined from 
measurements of x'( w ~~ > 0)), and the widths of spectral 
distributions of the resonant phonons Auj m k = = 
300 MHz according to the EPR data on the widths of 
the singlet-doublet transitions^, we found it necessary to 
introduce two different phonon lifetimes independent on 
temperature to describe the experimental data for the 
sample with the concentration c=0.104% (see FigsJ5]and 
[9]). Namely, we used T p h(w) = r g ~ 1/is for resonant 
phonons with frequencies corresponding to transitions 
between the electron-nuclear sub-levels of the ground elec- 
tronic doublet, and T p h(uS) = t s ~ 0.1/is for phonons 
with frequencies corresponding to transitions between 
the excited singlet and the ground doublet. Actually, 
the transitions between the electron-nuclear sublcvcls of 
the ground doublet arc narrower than the transitions be- 
tween the sublevels of the singlet and doublet crystal 
field states, and the lifetimes T g and t s may differ less 
than by order of magnitude. However, to fit the mea- 
sured frequency dependences of the susceptibilities in the 
samples with higher concentrations of holmium ions, we 
had to diminish the factors K m k by 1.48 (c=0.157%) and 
2.7 (c=0.27%) times as compared with these factors for 
c=0.104%. Thus, despite the increase of the concentra- 
tion of the paramagnetic ions, the phonon bottleneck ef- 
fect weakens due to broadening of the spectral bands of 
the resonant phonons and possible decrease of the phonon 
lifetimes. 

To illustrate the phonon bottleneck effect on the relax- 
ation rates, we present in Fig. IPTT a) the calculated rates 
in the sample with the concentration c=0.104% versus 
the magnetic field. These rates are condensed within the 
much narrower range as compared with the relaxation 
rates represented in Figs|3(a,b), the damping effect of 
the gaps at the Am = crossings is partly restored due 
to the strong suppression of the singlet-doublet transi- 
tion probabilities, but there are still no remarkable vari- 
ations of the lower branches at the crossing points which 
may be expected in accordance with the measured field 
dependences of the in- and out-of-phase susceptibilities 
represented in Figs 07] 

As it is seen in FigfT2]the magnetic field dependences 
of x' an d x" obtained without phonon bottleneck and 
cross-relaxation effects (with electron-phonon transition 
probabilities only, curves 1,2) have nothing in common 
with the experimental data (Figs(S][7|). When making use 
of electron-phonon transition probabilities renormalized 
by phonon bottleneck, we obtain the curves 3,4 show- 
ing some overlap with the measured ones (the specific 
humps and dips appear at the odd crossing points) but 
their "intensities" are lower than the experimental ones, 
and also there is no sign of corresponding features at 
even and half-integer crossings. However the results dra- 
matically change when the cross-relaxation is taken into 
account: in this case all crossing points manifest them- 
selves properly (the curves 5,6) and are very close to the 
experimental results. 

The cross-relaxation rates were calculated assuming 
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FIG. 12: The simulated ac-susceptibility of the LiYF 4 :Ho 3+ 
(0.27%) sample at the frequency 800 Hz and the tempera- 
ture 2 K, the magnetic field is declined from the c-axis by 
ldeg. a - the relaxation matrix contains only one-phonon 
transition probabilities calculated with the electron-phonon 
coupling constants presented in Table 1, b - the transition 
probabilities are renormalized due to the phonon bottleneck 
effect, c - the cross- relaxation terms are included. 



FIG. 11: The phonon bottleneck and cross relaxation effects 
on the relaxation rates, (a) - the relaxation rates at the tem- 
perature 2 K in the sample LiYF 4 :Ho 3+ (0.104%) calculated 
with the renormalized electron-phonon transition probabili- 
ties due to the bottleneck effect, (b) - results of calculations 
with the cross relaxation terms taken into account. 



the Gaussian line shape g CR — (27rA) _1 exp[— (uj pn — 
wim) /2A 2 ] of the spectral density of the energy reservoir 
corresponding to interactions between the holmium ions. 
It is enough to consider in (fT4"|) only three nonzero fit- 
ting parameters, k\2 and ^66, which determine the tran- 
sition probabilities within the manifold of the electron- 
nuclear sublevels of the ground doublet, and ^44, which 
determines rates of the singlet-doublet transitions. For 
the scaling factor S = 2.10 8 s -1 , the fixed values of 
ki2 = kae = 0.1 and £144 = 0.004 independent on the 
concentration of the holmium ions were used in all cal- 
culations. However, the dispersion A of frequencies of 
the cross-relaxation transitions depends on the concen- 
tration, it was estimated from the EPR lincwidths and 
corrected from a comparison of the simulated field de- 
pendences with the experimental data (A = 100, 120, 
140 MHz for the g%f line shape, and A = 185, 200, 240 
MHz for the g^ 2 R and g^P line shapes for samples with 
concentrations 0.104, 0.157 and 0.27%, respectively). As 
it is seen in FigfTTTb). the cross-relaxation rates play the 
dominant role at all crossing points. The most impor- 
tant result of the cross-relaxation processes is the appear- 
ance of the low frequency (10 2 — 10 3 s _1 ) branch in the 
spectrum of relaxation rates with the well resolved max- 



ima at the crossing points. It should be noted that the 
calculated field dependences of the relaxation rates may 
change remarkably when using another cross-relaxation 
line shapes (the Lorentz distributions, in particular)^, 
but to analyze this problem it is necessary to perform 
measurements in the samples with a more wide range of 
concentrations of paramagnetic ions. 

Small peaks at the half-integer crossings in the field de- 
pendences of x' (and the peaks and dips at these points 
in the field dependences of x") are induced by the cross- 
relaxation singlet-doublet transitions (they disappear in 
the simulated curves for £44=0). The value of the param- 
eter &44 used in simulations (the results are represented 
in FigsOZJ) coincides by an order of magnitude with the 
estimation based on the assumption of the prevailing role 
of the magnetic dipole-dipole interactions: 



fc 44 = 2^ca 6 (jif) 



here the sum is taken over all Y 3+ sites, a is the lat- 
tice constant, and c is the concentration of the impurity 
holmium ions. However, values of the parameters fci2 
and /egg mentioned above are much larger than the simi- 
lar estimations, and it is possible that the virtual phonon 
exchange contributes essentially into the interaction be- 
tween the holmium ions in the ground state, in particular, 
due to strong coupling with the dynamic lattice deforma- 
tions of B g symmetry. 

The peaks Cd in the high-field dependences of x' an d 
x" close to the magnetic field value B = 8AB in Fig[7] 
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correspond to the anti-crossing of the electron-nuclear 
sublevels of the excited singlet (see FigfJ}. These peaks 
appear due to a combined action of the random crys- 
tal field (described by the Hamiltonian Hs, see eq. (fT9")) ) 
and the cross-relaxation processes (singlet-doublet tran- 
sitions). So, we had a possibility to determine indepen- 
dently the &44 parameter from fitting intensities of Chi 
peaks, and parameters of the Hamiltonian Hs from fit- 
ting intensities of the Cd peaks. The best fit was achieved 
when the parameters B 2 , B^ 2 from Ref[^ (see section III 
above) were increased by a common factor of 1.5, 1.8 
and 3.3 for samples studied in this work with holmium 
concentrations 0.104%, 0.157% and 0.27%, respectively, 
and natural abundances of Li isotopes. It should be 
noted that we have not succeeded in exact description of 
measured real and imaginary parts of the dynamic sus- 
ceptibilities at any frequency, temperature and external 
magnetic field values. In particular, we obtained correct 
positions of the maxima in the frequency dependences 
of the out-of phase susceptibilities, but the calculated 
maximum values are higher than the measured data (see 
Fig. 9). The discrepancies between theory and experi- 
ment in Fig. 7 are the consequences of this lack of the 
theory. However, the specific behavior of the suscepti- 
bility in the vicinity of high field anti-crossings is repro- 
duced by calculations as well. To derive a more elab- 
orated model, we need additional experimental data at 
higher frequencies. 

It follows from the simulations that an additional nar- 
row peak should split out from the peak Cd in the plot 
of x' versus the magnetic field (see the curve 5 in Fig fT2^) 
in the transverse magnetic field, in particular, due to a 
rather small misalignment of the sample. Because such 
a peak was not observed, we had good reasons to believe 
that the samples were oriented with the accuracy better 
than 1 degree. 

Using the model parameters determined from the anal- 
ysis of the frequency and field dependences of x' an d x" \ 
we obtained a satisfactory description of the temperature 
dependences as well (see Fig flO|) . The transformations of 
the magnetic field dependences of the out-of-phase sus- 
ceptibility at the fixed frequency in the vicinity of the 
crossing points with temperature fFigs l5l7|) or concen- 
tration (FigdJJ) can be explained as a result of strong 
variation of relaxation rates with the magnetic field at 
these points. If the effective relaxation rate exceeds the 
frequency of the ac field, x" has a dip, and if the re- 
laxation rate is less than the frequency, x" has a peak 
at the crossing point. The corresponding "peak to dip" 
transformations can be easily recognized in Figs EH 

VII. CONCLUSION 

We conclude that the microscopic model of a linear 
non-resonant response of the electron-nuclear subsystem 



in the highly diluted paramagnetic crystal LiYF4:Ho 3+ 
on the weak ac field has been derived. The model op- 
erates with a few phenomenological parameters intro- 
duced to account for a finite heat capacity of the res- 
onant phonons and cross-relaxation processes. The inde- 
pendent on temperature finite phonon lifetimes of about 
1 /us, comparable with the phonon time of flight between 
the sample boundaries, have been determined from fit- 
ting the results of calculations to the measured frequency 
dependences of out-of-phase susceptibilities in the sam- 
ples with different holmium concentrations. The parame- 
ters of the interaction between paramagnetic ions, which 
determine cross-relaxation rates, were obtained from fit- 
ting the simulated dependences of in-phase and out-of- 
phase susceptibilities on the external magnetic field. The 
obtained value of the parameter k44, which defines co- 
tunneling processes at half-integer crossing points in the 
spectrum of a pair of Ho 3+ ions, agrees with the estima- 
tion based on the dipole-dipole mechanism of inter-ion 
coupling. However, to present definite conclusions about 
the most important terms in the Hamiltonian of inter- 
action between the holmium ions, which are responsible 
for the cross-relaxation processes at integer crossing and 
anti-crossing points, it is necessary to derive a more elab- 
orated theory of cross-relaxation than the high tempera- 
ture approach of Ref. [5] . 

The model containing single-ion crystal-field includ- 
ing random strains, electron-phonon transitions includ- 
ing bottleneck, hypcrfine and cross-spin interactions, is 
comprehensive and can be used to predict the detailed 
magnetization dynamics in LiYF4:Ho 3+ crystals at ul- 
tra low temperatures and at higher frequencies. Note 
that half-integer transitions were observed at elevated 
temperatures (above 1.5 K) only in ac susceptibility ex- 
periments or at fast sweeping field in the low tempera- 
ture (below 0.1 K) magnetization measurements where 
the bottleneck effects stabilized an effective spin-phonon 
temperature of half a Kelvin . In the present the- 

oretical approach half-integer transitions result from res- 
onant cross-relaxation transitions involving sublevels of 
the lower singlet with the activation energy of ~10 K, 
which is satisfactory to fit the experimental data taken 
at liquid helium temperatures. A more detailed analysis 
of the magnetization dynamics in the sweeping fields will 
be presented in a separate paper. The theory may be ad- 
ditionally tested by measurements of the ac susceptibility 
and NMR in a tilted magnetic field. Finally, this theo- 
retical approach may be also expanded to other magnetic 
systems, in particular, the single molecule magnets. 
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